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INTRODUCTION 


In  simplest  terms  EM-TRANAIR  is  a  computer  program  for  the  solution  of 
Maxwell’s  equations  in  three  dimensions.  The  ability  to  execute  the  code 
for  the  Aerodynamic  Methods  Group  marks  the  beginning  of  the  investigation 
into  low  observable  technology.  EM-TRANAIR  was  developed  by  the  Boeing 
Military  Airplane  Company  and  delivered  to  the  USAF  in  May  1987. 


"EM-TRANAIR  uses  a  variational  principle  to  characterize  Maxwell's 
field  equations  and  specified  boundary  conditions  as  an  extremum  problem  in 
the  calculus  of  variations.  A  multiply-preconditioned  Krylov  space  method 
called  GMRES  (generalized  minimum  residual  algorithm) ,  with  both  an  extreme¬ 
ly  powerful  sparse  solver  preconditioner  and  an  exterior  Helmholtz  solver 


preconditioner,  provides  convergence  to  machine  error  in  ten  to  twenty 
iterations.  " 


Some  schemes  used  to  calculate  low  observable  parameters  initially 
determine  the  point  sources  on  the  surface  grid  generated  by  the  electric  or 
magnetic  field.  From  these  point  sources  the  far  field  parameters  are 
calculated  with  a  numerical  analysis.  Through  a  different  approach  using 
prescribed  mathematical  tools,  EM-TRANAIR  uses  sources  that  are  fundamental 
unknowns  for  the  constant  coefficient  differential  operator  governing  the 
far  field,  rather  than  the  electric  or  magnetic  field  at  grid  points.  The 
solution  is  peculiar  to  the  numerical  analysis  of  Maxwell's  equations. 

Thi6  report  is  in;ended  to  relate  information  concerning  manipulation 
and  execution  of  the  EM-TRANAIR  code.  The  process  of  code  installation 
became  as  involved  as  the  actual  theory  and  mathematical  analysis  employed 
to  achieve  solution.  Hopefully  the  following  content  will  prove  to  be  a 
useful  tool  to  troubleshoot  the  EM-TRANAIR  software. 

EM-TRANAIR  consists  of  two  main  programs  supported  by  five  libraries. 
The  two  main  programs  are  an  input  processor,  V1INP,  and  a  solver/ 
post-processor,  V1S0L.  The  input  processor  is  extremely  extensive  in  length 
which  makes  it  difficult  to  troubleshoot  for  errors.  It's  quite  time 
intensive  when  it  comes  to  gaining  familiarity. 

The  library,  input  and  solver  files  are  first  saved  as  program  listings 
and  binary  files  using  UPDATE.  "UPDATE  is  a  line  oriented  test  editor  for 
maintaining  programs  in  the  form  of  source  code,  as  well  as  other  types  of 
text  data.  UPDATE  creates  and  modifies  program  libraries  and  produces 
output  that  can.be  used  as  input  to  other  programs,  particularly  compilers 
and  assemblers.  " 

The  library  files  are  needed  so  that  pieces  may  be  pulled  from  them 
when  the  source  decks  for  the  input  processor  and  the  solver  are  created 
using  SEGLDR.  The  reason  for  this  methodology  is  to  economize  the  available 
central  memory  during  code  execution.  Use  of  SEGLDR  results  in  the  creation 
of  a  nonsegmented  executable  program.  One  for  the  input  processor  called 
EMINP  and  one  for  the  solver  called  EMSOL  will  be  created.  Nonsegmented 
programs  are  those  that  have  all  of  their  code  continually  memory  resident. 
"SEGLDR  is  an  automatic  loader  for  code  produced  by  language  processors.  " 
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The  code  delivered  by  Boeing  was  compiled  using  CFT  FORTRAN  compiler  on 
a  CRAY  computer  with  8  million  words  of  available  central  memory  under  the 
COS  operating  system.  These  three  characteristics  (chosen  compiler,  size  of 
central  memory  and  chosen  operating  system)  are  key  factors  influencing 
successful  code  execution. 

Execution  problems  needed  to  be  resolved  to  utilize  EM-TRANA1R  and 
obtain  useful  numerical  RCS  calculations.  Compiler  problems  such  as  the 
ability  to  pass  argument  values  between  source  programs  and  subroutines  are 
addressed.  Hardware  limitation  such  as  the  size  of  the  central  core  memory 
is  an  important  factor  in  what  configurations  would  result  in  converged 
solutions.  Manipulations  within  the  FORTRAN  codes  such  as  setting  initial 
values  for  variables  and  memory  allocation  for  vector  arrays  are  also 
discussed.  Some  familiarity  with  successful  numerical  RCS  calculations 
resulted  in  the  presentation  of  a  few  "rules  of  thumb."  These  and  other 
installation  and  execution  problems  are  presented  to  help  the  reader  with 
EM-TRANAIR. 

This  document  is  not  a  dissertation  on  theory  and  code  development  but  a 
helpful  primer  to  EM-TRANAIR  installation  at  WRDC  and  a  guide  to  trouble¬ 
shooting  confronted  execution  and  processing  problems. 


2 


INPUT  DATA/ SOURCE  DECK 


GFLIB,  EMLIB,  GPLIB,  TRLIB  and  EXSOL  are  the  libraries  that  support  the 
input  processor  and  the  solver.  GFLIB  calculates  the  Green's  function. 

EMLIB  contains  the  EM  (electromagnetic)  subroutines.  GPLIB  contains  the 
general  purpose  subroutines.  TRLIB  contains  supplemental  subroutines. 

EXSOL  contains  the  exterior  Helmholtz  solver.  The  subroutines  within  the 
libraries,  the  input  processor  and  the  solver  are  contained  in  Appendix  A. 

Appendix  B  shows  a  listing  generated  when  the  input  processor  V1INP  was 
initially  compiled  using  CFT77.  Because  of  the  difference  between  the  CFT 
and  CFT77  compilers  a  number  of  bugs  were  created  and  needed  to  be  resolved 
for  proper  compilation.  Equivalence  and  common  statements  left  a  number  of 
errors  as  can  be  seen  in  the  listing  in  Appendix  B.  The  CFT77  compiler 
didn't  like  the  way  variables  were  being  passed  from  the  source  code  back 
and  forth  to  the  subroutines.  The  call  and  subroutine  statements  involving 
these  passed  variables  were  extended  to  fully  cover  all  the  arguments.  In 
this  way  the  equivalence  and  common  statements  could  be  deleted  from  the 
code  in  certain  places.  The  best  and  recommended  way  to  find  the  changes 
would  be  to  compare  the  original  deck  with  the  executable  deck.  It  may  be  a 
time  consuming  task,  but  the  decks  are  just  too  extensive  to  create  listings 
and  the  number  of  changes  were  numerous.  The  listing  such  as  in  Appendix  B 
will  provide  the  information  so  one  can  zero  in  on  an  area  of  interest. 
Remember,  the  listing  in  Appendix  B  can  be  generated  for  any  code  with  ON=H 
in  the  CFT77  job  control  statement. 

Another  problem  encountered  with  compiling  on  CFT77  was  the 
dimensioning  of  arrays.  The  pointers  for  the  arrays  weren't  being 
calculated  as  integers  but  real  numbers .  These  problems  were  fixed  by 
insuring  arithemic  operations  dealing  with  array  pointers  resulted  in 
integers. 

The  final  problem  worth  mentioning  dealt  will  boolean  algebra 
operations.  Only  the  CFT77  can  do  boolean  manipulations.  Makes  you  wonder 
how  the  Boeing  folks  got  around  this  one  using  CFT?  The  bottomline  is  that 
an  executable  deck  containing  boolean  algebra  requests  can  only  be  compiled 
using  CFT77  to  be  successful  at  WRDC. 

UNICOS  operating  system  requires  the  use  of  the  CFT77  compiler.  Half 
the  battle  is  won  switching  from  COS  to  UNICOS  because  hopefully  most  of  the 
bugs  have  been  resolved.  The  code  has  been  made  CFT77  compile  successful 
although  being  originally  compiled  on  CFT. 

Just  a  little  hint  on  troubleshooting,  the  $DEBUG  and  dump  files  are  a 
great  aid  in  getting  traceback  maps.  These  maps  help  in  following  paths  to 
that  place  in  the  code  that  tripped  the  execution  error.  Using  COS  or 
UNICOS,  various  flags  are  turned  on  in  the  job  control  language  statements 
to  establish  symbol  files  and  maps  for  $DEBUG. 

Happy  hunting. 


3 


V1INP  INPUT  PROCESSOR 


The  input  processor  is  the  heart  of  EM-TRANAIR.  The  code  contains 
approximately  27,000  lines  making  path  following  to  debug  problems  quite  a 
task.  The  nonsegmented  structure  helps  in  your  searches  for  execution 
problems  because  the  numerous  subroutines  plus  the  libraries  subroutines  act 
as  blaze  markings  along  the  pathway.  The  imbedded  error  flags  and  output 
messages  are  also  another  big  plus. 

The  processor  reads  in  the  case  test  configuration  and  does  extensive 
presolving  tests  and  comparisons  to  avoid  glitches  when  it's  time  to  use  the 
solver.  "The  V1INP  program  reads  a  description  of  a  configuration  and  some 
incidence  conditions  describing  a  plane  wave  and  computes  coefficients  of  a 
discrete  operator  on  a  Cartesian  grid  which  approximates  Maxwell's  equations 
and  appropriate  configuration  boundary  conditions. 

The  functional  operations  of  the  input  processor  are: 
o  Process  User  Input  Data 

o  Assemble  Configuration  and  Test  for  Validity 
o  Define  Coupling  with  Cartesian  Grid 
o  Define  Locations  of  Unknowns 
o  Compute  Discrete  Operators 
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o  Communicate  Problem  Definition  to  Solver,  V1S0L  " 

The  first  problem  and  probably  the  most  recurring  is  central  memory 
management.  Because  the  size  available  was  limited  to  under  2  million  words 
of  memory  the  choice  of  a  test  case  was  highly  restrictive.  Various  cases 
run  by  Boeing  are  discussed  in  reference  [1]  and  one  of  these  would  be  a 
good  example  to  compare  against  during  code  installation  execution  runs.  It 
became  fairly  apparent  that  only  the  sphere  would  become  the  likely 
candidate. 

The  key  driver  is  choosing  a  size  for  the  parameter  NDUMSS.  The 
default  values  or  really  the  values  set  in  the  delivered  code  established 
for  NDUMSS  and  other  memory  management  parameters  are  listed  in  Table  1. 
NDUMSS  sets  the  amount  of  scratch  storage  needed  for  the  input  processor  to 
do  its  compilations  and  manipulations.  Using  NDUMSS  set  to  1.1  million 
words  really  limits  the  amount  of  remaining  central  memory  for  the  other 
parameters  when  less  than  2  million  words  total  is  available.  With  this  in 
mind  and  discussions  with  Boeing,  NDUMSS  was  set  to  800,000  words  and  the 
sphere  would  be  the  only  test  case  that  might  fit. 

The  sphere  configuration  was  still  to  large  to  fit  on  the  CRAY  using 
COS.  The  solution  was  to  pull  out  every  other  point  of  the  surface  grid. 
This  resulted  in  an  executable  test  configuration.  Using  COS  and  less  than 
2  million  words  central  memory  resulted  In  a  really  simple  geometry  for  the 
initial  code  installation.  Memory  management  is  highly  dependent  upon  input 
configuration. 


4 


Variable 

Value 

Description 

MXNETT 

225 

Maximum  number  of  networks 
in  configuration 

MAXPTS 

7500 

Maximum  number  of  mesh  points. 

NGRD 

140481 

Maximum  number  of  grid  points. 

MAXSRC 

14048 

Maximum  number  of  sources, 
any  one  component 

NDUMSS 

1.1M 

Total  scratch  storage 

MAXSRCH 

100 

Maximum  number  of  search  directions 
(GMRES)  any  one  component 

MAXMAT 

100 

Maximum  number  of  different  materials 

MXLPAN 

150 

Maximum  number  of  panels  in  a  grid  box 

MXNTPN 

2000 

Maximum  number  of  panels  in  one  network 

MAXNB 

14048 

Maximum  number  of  B  points  (not  used) 

MXIABT 

3200 

Maximum  of  pairwise  abutments 

MXNABT 

1000 

Maximum  number  of  abutments 

MXEIAB 

20 

Maximum  number  of  edges  in  an  abutment 

MXFDSG 

1200 

Maximum  number  of  fundamental  segments 

MXEMPT 

6500 

Maximum  number  of  edge  mesh  points 

MXEDMP 

400 

Maximum  number  of  edge  mesh 
points  per  network 

MXNAI 

350 

Maximum  number  of  abutment  intersections 

MXNPEC 

100 

Maximum  number  of  pts  in  an  eq.  class 

MNODMX 

100 

Maximum  number  of  nodes 
at  abutment  intersection 

NSEGMX 

60 

Maximum  number  of  segments 
at  abutment  intersection 

MAXKMS 

19500 

Maximum  number  of  micro-edge 
abutment  descriptions 

TABLE  1:  EM-TRANAIR  Parameter  Declarations 
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Another  area  that  dictated  concern  dealt  with  the  initialization  of 
code  parameters  and  variables.  For  some  unknown  reason  the  original  decks 
for  the  input  processor  and  the  solver  wouldn't  execute  properly  with  the 
initial  value  sets  for  the  variables.  For  the  processor  and  solver,  source 
program  AAINPUT  and  AA3DS  respectively,  would  call  up  a  subroutine 
containing  the  initial  value  sets.  Appendix  C  contains  the  listing  of  the 
initial  value  set  for  the  input  processor. 

To  get  around  this  execution  problem  with  the  passing  of  the  initial 
values  the  contents  of  the  initial  values  subroutine  was  imbedded  into  the 
source  program.  This  way  the  transfer  with  the  call  and  return  was 
eliminated.  If  the  BLOCK  DATA  statement  was  left  in  the  source  program  the 
execution  would  still  stop.  So  the  BLOCK  DATA  statement  was  also 
eliminated. 

When  the  "new"  deck  with  the  imbedded  initialization  values  is  run  a 
caution  is  returned  during  the  CFT77  compile  operation.  The  caution  is  to 
reveal  that  with  the  initialized  values  imbedded  in  the  source  program 
there's  no  guarantee  that  the  initialized  value  will  be  kept  throughout  the 
program.  In  the  original  deck  the  initialize  subroutine  is  only  called  once 
so  its  pretty  safe  to  say  no  problems  are  encountered  due  to  initialization 
imbedded  in  the  program  source  deck.  The  imbedded  initialization  to  the 
input  processor  V1INP  program  source  deck  AAINPUT  is  also  shown  in  Appendix 
C. 


As  stated  previously,  central  core  memory  management  was  a  recurring 
problem.  After  getting  over  the  variable  initialization  problem  the  input 
processor  still  wouldn't  execute  too  far  down  into  the  deck.  A  bit  of  time 
was  consumed  troubleshooting  and  the  best  means  of  solution  was  to  track  the 
values  of  various  variables  using  write  statements  to  print  out  values  on 
FT06.  It  was  discovered  that  the  surface  grid  network  for  the  conducting 
sphere  had  twice  as  many  points  as  there  was  space  provided  for  in  the 
vector  arrays.  There  weren't  any  imbedded  helpful  error  messages  in  V1INP 
to  point  at  a  fix  to  the  problem.  Since  the  size  of  central  memory  was  set, 
the  number  of  points  on  the  surface  grid  had  to  be  decreased.  Every  other 
point  was  eliminated  on  the  surface  grid  input  data.  The  original  surface 
network  grid  and  final  network  grid  for  the  conducting  sphere  configuration 
are  shown  in  Appendix  D.  Using  the  UNICOS  operating  system  Instead  of  COS 
because  of  more  available  central  memory  should  alleviate  any  problems  with 
surface  grid  sizes. 

One  oversight  to  avoid  when  trying  to  follow  tracebacks  concerns  calls 
to  subroutines  that  won't  be  found  in  the  listings  of  Appendix  A.  It's  not 
that  they  don't  exist  but  they  are  CRAY  libraries.  A  couple  of  discovered 
examples  where  CRAY  subroutines  are  called  upon  are  as  follows. 

The  library  deck  GPLIB  which  contains  the  general  purpose 
subroutines  also  has  a  call  statement  for  WOPEN.  When  the  subroutine 
listings  for  the  libraries  in  Appendix  A  are  reviewed  it  will  be  discovered 
that  there  is  no  listing  for  WOPEN.  This  is  an  internal  CRAY  library 
subroutine  that  opens  a  word-addressable,  random  access  dataset. 
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The  same  can  be  discovered  for  the  routine  SCOPY  found  in  the 
input  processor  deck,  subroutine  ABTIDN.  SCOPY,  again  an  internal  CRAY 
library  subroutine,  is  used  to  copy  a  real  or  complex  vector  into  another. 

Tracebacks  greatly  aid  in  path  following  troubleshooting.  If  the 
path  appears  to  call  on  an  unlisted  subroutine  not  found  in  Appendix  A,  it's 
probably  a  CRAY  internal  library  subroutine. 

The  various  problems  encountered  by  the  author  to  get  the  input 
processor  V1INP  to  execute  have  been  discussed.  If  any  hints  or  insights 
were  revealed  to  the  next  individual  working  with  EM-TRANAIR  the  reward  was 
well  worth  the  effort.  Use  this  information  to  work  the  changes  made  to  the 
original  decks  if  different  test  configurations  result  in  unexpected 
execution  errors  and  unrealistic  computational  results. 


7 


VI SOL  SOLVER 


With  the  input  processor  being  the  heart  of  EM-TRANAIR,  the  solver, 
V1SOL  is  the  life  support.  Extensive  preprocessing  with  the  input  processor 
establishes  an  input  configuration  that  predictably  runs  through  the  solver 
error  free.  By  comparison  the  input  processor  contains  approximately  27,000 
lines  of  code  while  the  solver  contains  approximately  3,000  lines  of  code. 


The  solver,  V1S0L,  "relies  heavily  on  two  special  purpose  libraries  to 
provide  the  discrete  Green's  function  and  the  exterior  Helmholtz  Solver 
(GFLIB  and  EXLIB,  respectively).4-1" 


"The  V1SOL  program  iteratively  solves  the  set  of  discrete  equations 
which  approximate  Maxwell's  equations  and  describes  the  scattering  of 
electromagnetic  energy  from  some  material  object.  The  solution  process  can 
be  divided  into  the  tasks  of: 


o  Initialization 
o  Data  Entry 
o  Overhead  Operations 

r  4 1 

o  The  Iterative  Solution  Process  " 

As  with  the  input  processor,  the  problem  of  initialization  of 
parameters  was  also  encountered  with  the  solver.  Two  subroutines,  INIT  and 
INICOM  (shown  in  Appendix  E) ,  In  the  original  V1S0L  deck  were  used  to  pass 
parameter  initial  values  into  the  source  program  AA3DS.  The  source  program 
AA3DS  would  call  on  INIT  which  would  call  on  INICOM.  As  expected,  the 
original  deck  wouldn't  execute.  To  resolve  the  problem,  subroutine  INICOM 
was  imbedded  into  subroutine  INIT  which  in  turn  was  imbedded  into  the  source 
program  AA3DS. 


The  results  are  shown  In  the  second  listing  of  Appendix  E.  This 
listing  is  the  revised  section  of  the  source  program  AA3DS  that  contains  the 
imbedded  parameter  initialization. 

The  only  other  execution  error  encountered  with  the  V1S0L  was  due  to  a 
couple  of  vector  arrays  being  erased  from  allocated  central  memory.  The 
vector  array  was  then  needed  for  manipulation  further  into  the  code  but  was 
no  longer  available.  This  problem  occurs  in  subroutines  CLUSS  and  CLUSM 
which  are  found  in  V1S0L  and  called  from  the  source  program  AA3DS.  The 
subroutines  are  listed  in  Appendix  F  along  with  a  "dummy"  subroutine  called 
IEQUAL. 


Vector  arrays  C3  and  KCB  are  the  ones  that  are  inadvertently  erased. 

In  subroutine  CLUSS,  vector  array  KCB  is  equated  with  vector  array  NC.  When 
this  is  done  vector  array  KCB  becomes  unlisted  on  central  memory  because 
that  allocation  of  memory  space  has  become  known  as  NC.  This  is  where  the 
IEQUAL  "dummy"  subroutine  comes  In.  By  passing  the  vector  arrays  CB  and  KCB 
as  arguments  their  memory  allocation  remains  intact.  If  this  isn't  done  the 
computer  can't  find  the  arrays  because  essentially  they  don't  exist.  But 


8 


with  IEQUAL,  vector  arrays  CB  and  KCB  still  retain  their  values  and  memory 
allocations. 


Why  are  both  vector  arrays,  CB  and  KCB,  passed  as  arguments  to 
subroutine  IEQUAL?  Depending  on  subroutine  CLUSS  and  CLUSM,  there  are 
statements  within  these  subroutines  that  equate  CB  and  KCB  to  other  vector 
arrays.  Note  Appendix  F  is  what  code  exists  in  the  revised  V1S0L  deck. 
IEQUAL  didn't  appear  in  the  original  V1S0L  deck. 

Installing  the  solver  was  a  small  hurdle  compared  to  installing  the 
input  processor.  Hopefully,  the  insights  given  to  troubleshooting  these 
execution  errors  will  aid  in  any  future  troubleshooting.  These  solutions 
may  be  a  place  to  start  an  investigation  for  future  solver  software  problems 
if  encountered. 
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CONFIGURATION  EXECUTION  PROCESSING 


Before  the  SSD  on  the  CRAY  was  committed  as  an  external  resource  for 
UNICOS  a  number  of  configurations  where  submitted  for  execution.  This 
section  deals  with  a  few  "rules-of-thumb"  that  were  discovered  with  this 
opportunity  to  test  the  code  while  the  COS  operating  system  was  still 
available.  These  rules-of-thumb  are  still  applicable  with  the  UNICOS 
operating  system  because  they  are  peculiar  to  EM-TRANAIR. 

When  execution  starts  in  the  solver,  V1S0L,  various  parameters  are 
computed  and  output  to  FTC6.  Two  of  these,  delta  (for  X,  Y  and  Z)  and  the 
number  of  mesh  points  per  wavelength,  lambda,  are  especially  important. 
Delta  should  be  within  the  range 

.2  <  delta  <  5 

and  the  number  of  mesh  points  per  wavelength,  lambda,  should  be  within  the 
range 


5  <  poiuts/A<  10 

If  you  stay  within  these  limits  the  numerical  results  are  reliable. 

The  Green's  function  is  calculated  within  V1S0L.  Information 
concerning  the  Green's  function  is  output  to  F106.  Because  of  some  bug  in 
the  cede  a  flag  stays  on  that  causes  error  tolerances  to  be  printed  even  if 
the  calculations  are  within  accurate  solution  boundaries.  Note  that  if  the 
tolerances  output  to  FT06  are  zeros;  rest  assured  the  Green’s  function  is 
accurate. 

Beware  of  thin  geometries  such  as  disks  and  flat  plates.  The  delta 
aspect  ratio  in  the  "thin"  direction  may  lead  to  problems  in  getting  a 
converged  solution  to  Maxwell's  equations. 

A  "thin"  configuration  resulting  in  a  solution  is  greatly  influenced  by 
the  chosen  frequency  used  for  the  incidence  radiation.  If  a  low  frequency 
is  specified,  you  won't  have  many  mesh  points  per  wavelength,  lambda,  in  the 
"thin"  direction.  If  a  high  frequency  is  specified  you’ll  be  playing  the 
aspect  ratio  deltas  against  the  calculated  discrete  wave  number.  A  possible 
solution  to  "thin"  geometries  may  be  to  scale  up  the  geometry  ($REFerence 
length  parameter)  so  the  "thin"  dimension  is  no  longer  a  problem. 

Finally,  be  careful  with  the  $B0X  definition  parameters.  Be  sure  the 
number  of  mesh  points  in  the  computational  grid  is  of  the  type 

(2a)(3b)(5C)+l 

Where  a,  b  and  c  are  integers  and  the  resulting  number  is  odd. 

These  "rules-of-thumb"  will  be  helpful  in  achieving  converged 
solutions.  Experience  with  the  code  is  the  best  medicine  but  the 
rules-of-thumb  are  a  good  heads  up  before  diving  in. 
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POST  RCS  RESULTS 


Before  switching  the  SSD  over  to  UNICOS  a  few  weeks  were  given  where  4 
million  words  of  central  memory  were  available  with  COS.  This  was  an 
excellent  opportunity  to  try  some  EM-TRANAIR  runs  before  having  to  change  to 
UNICOS.  This  was  the  chance  to  gain  some  familiarity  with  the  code  and 
become  better  assured  things  were  executing  properly. 

Various  configurations  were  tried  and  some  were  successful  and  others 
ran  into  problems  as  were  outlined  in  che  previous  section  concerned  with 
configuration  execution  processing.  Presented  here  are  some  actual  post  RCS 
graphics  that  compare  results  to  those  conducted  by  Boeing  and  an  actual 
comparison  to  an  experimental  run  taken  at  one  of  ASD/WRDC  radar  ranges. 
These  preliminary  results  are  encouraging  and  reveal  the  acceptable 
reliability  concerning  the  software  code  EM-TRANAIR. 

The  first  configuration  used  was  a  conducting  sphere.  This 
configuration  is  the  same  as  used  by  Boeing  for  code  validation  and  a 
precise  definition  is  presented  in  reference  [1,4].  Figure  1  was  taken  from 
reference  [1].  Figure  2  shows  the  results  for  the  same  conditions  as  Figure 
1  but  being  an  execution  run  on  WTAFB's  CRAY  X-MP/216.  A  total  RCS 
measurement  for  vertical  polarization  incidence  at  2  degree  interval 
bistatic  stations  is  calculated.  A  fifth  order  polynomial  curve  fit  was 
used  to  create  the  smooth  fit  through  90  bistatic  stations.  The  agreement 
is  excellent. 

The  frequency  of  the  incident  radiation  is  the  driver  that  regiments 
acceptable  and  non-acceptable  configuration  runs  as  far  as  number  of  mesh 
grid  points  and  points  per  wavelength  lambda.  Numerous  configurations  were 
tried  to  get  a  better  feel  for  the  driving  parameters  and  the  number  of 
failures  exceeded  the  number  of  successes.  The  size  of  an  acceptable 
geometry  is  determined  by  the  incidence  radiation  frequency.  This 
relationship  is  the  key  to  an  accurate  solution  of  Maxwell's  equations  using 
Green's  function. 

Initial  runs  were  executed  at  high  incident  frequencies  from  9  to  16 
GHz.  The  geometries  being  used  just  weren't  compatible  to  the  desired 
incidence  frequencies.  Possibly,  the  problem  was  being  to  close  and  in  some 
cases  being  involved  within  the  resonance  region.  Looking  over  the  evidence 
mostly  the  Rayleigh  region  appears  to  lead  to  code  execution  success. 

The  bottom  line  for  the  preliminary  lessons  learned  revealed  that 
staying  within  the  Rayleigh  region  brought  about  accurate  solutions  and  this 
is  the  norm  for  the  configurations  used  by  Boeing  for  validation.  So  more 
configurations  were  executed  to  try  to  mimic  Boeing's  results. 

Moving  up  in  incidence  frequency.  Figure  3,  taken  from  reference  [1]  is 
Boeing's  run  with  a  conducting  sphere  at  a  medium  frequency  as  determined  by 
the  diameter  of  the  sphere.  The  execution  results  on  the  CRAY  X-MP/216  for 
the  same  configuration  used  for  Figure  3  are  presented  in  Figure  4.  A  total 
RCS  measurement  for  vertical  polarization  incidence  at  2  degree  interval 
bistatic  stations  was  calculated.  A  fifth  order  polynomial  curve  fit  was 
used  to  create  the  smooth  fit  through  90  bistatic  stations.  The  agreement 


11 


O'n  o'?i  col  0'6  0'9  o->  c:  c:  ct-  o-v- 

i wS3C 1  33a 


(wsaa)  sou 


12 


FIGURE  1:  Vertical  Polarization  Bistatic  FIGURE  2:  Vertical  Polarization  Bistatic 

RCS  of  Low  Frequency  Sphere  RCS  of  Low  Frequency  Sphere 

(Boeing)  (CRAY  X-MP/216) 


looks  good  with  the  possible  exceptions  of  the  deep  troughs.  More  bistatic 
calculations,  say  at  every  degree  interval,  would  produce  a  better  curve. 

Having  moved  up  in  incidence  frequency,  a  look  was  taken  to  see  what 
type  of  agreement  was  found  at  the  same  incidence  frequency  but  with 
horizontal  polarization  of  the  incidence  electric  field.  Figure  5  was  taken 
from  reference  [1].  Figure  6  shows  the  results  for  the  same  conditions  as 
Figure  5  but  being  an  execution  run  on  the  CRAY  X-MP/216.  A  total  RCS 
measurement  for  horizontal  polarization  incidence  at  2  degree  interval 
bistatic  stations  was  calculated.  A  fifth  order  polynomial  curve  fit  was 
used  to  create  the  smooth  fit  through  90  bistatic  stations.  The  agreement 
is  excellent. 

Finishing  up  the  conducting  sphere  configuration  successful  results,  is 
a  look  at  monostatic  RCS  measurements.  Figure  7  shows  the  monostatic  RCS 
measurement  for  a  6-inch  diameter  conducting  sphere  being  radiated  by  a 
vertical  polarized  electric  field  at  10  GHz.  This  actual  experimental  run 
is  presented  to  show  how  the  total  RCS  magnitude  is  a  constant  over  the 
surface  of  the  conducting  sphere.  A  good  numerical  result  using  EM-TRANAIR 
to  compare  to  the  experimental  measurements  in  Figure  7  was  not  achieved.  A 
good  guess  as  to  the  inaccuracies  in  the  Green's  function  solution  of 
Maxwell's  equations  may  be  due  to  the  high  frequency  incidence  radiation  and 
the  small  diameter  geometry.  A  converged  solution  couldn't  be  obtained. 

The  sphere  configuration  used  as  Boeing's  validation  case  was  used  to  obtain 
a  numerical  monostatic  RCS  solution  (the  low  frequency  case  with  ka-3.0). 

The  important  comparison  between  the  numerical  results  and  the  experimental 
results  in  Figure  7  is  the  constant  magnitude  of  RCS  over  the  sphere 
surface.  Figure  8  shows  the  total  RCS  measurement  for  vertical  polarization 
incidence  at  2  degree  interval  monostatic  stations.  A  fifth  order 
polynomial  curve  fit  was  used  to  create  a  smooth  fit  through  45  monostatic 
stations.  Figure  9  shows  the  total  RCS  measurement  for  horizontal 
polarization  incidence  at  2  degree  interval  monostatic  stations.  Once 
again,  a  fifth  order  polynomial  curve  fit  was  used  to  create  a  smooth  fit 
through  45  monostatic  stations.  The  monostatic  RCS  numerical  results  with 
EM-TRANAIR  definitely  possess  a  constant  RCS  magnitude  over  the  conducting 
sphere  surface  for  both  a  horizontal  and  vertical  polarized  incidence 
electric  field. 

The  last  configuration  and  numerical  RCS  result  to  be  discussed  was  an 
attempt  to  get  RCS  information  on  an  8-inch  diameter,  .25-inch  thick, 
conducting  circular  disk.  The  reason  this  configuration  was  chosen  is  due 
to  the  experimental  results  made  available  from  one  of  WPAFB  radar  ranges. 
These  experimental  results  using  an  Incidence  radiation  (vertical 
polarization)  at  6  GHz  are  shown  in  Figure  10.  The  EM-TRANAIR  numerical 
results  are  shown  in  Figure  11  for  the  conducting  circular  disk.  A  total 
RCS  measurement  for  vertical  polarization  incidence  at  2  degree  interval 
monostatic  stations  are  calculated.  A  fifth  order  polynominal  curve  fit  was 
used  to  create  a  smooth  fit  through  45  monostatic  stations.  The  Radar  Range 
(Figure  10)  angle  in  degrees  from  90  to  180  degrees  should  be  compared  to 
the  X-MP/216  (Figure  11)  angle  in  degrees  from  0  to  90  degrees.  This 
comparable  range  runs  from  a  face  on  monostatic  view  of  the  disk  revolving 
to  an  edge  on  monostatic  view.  Magnitude  in  DB  from  Figure  10  is  the  same 
as  Total  RCS  in  Figure  11.  The  comparison  between  the  experimental  and 
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FIGURE  7:  ASD/WRDC  Radar  Range  6- inch  Diameter  Conducting  Sphere  Test  Case 
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FIGURE  8:  Vertical  Polarization  Monostatic  FIGURE  9:  Horizontal  Polarization  Monostatic 

RCS  of  Low  Frequency  Sphere  RCS  of  Low  Frequency  Sphere 

(CRAY  X-MP/216)  (CRAY  X-KP/216) 
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FIGURE  10:  ASD/WRDC  Radar  Range  8-inch  Diameter 
25-inch  Thick  Circular  Disk  Test  Case 
Vertical  Polarization  Monostatic  RCS 
at  6  GHz 

015k  vtPT  POL  6  GHz 


I  15.0  2G.0  L5.0  3-j.O  35.0  1C.0  15. 0  50.0  55.0  bC.u  $5.G 
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FIGURE  II:  8-inch  Diameter  .25-inch  Thick 

Circular  Disk  Vertical  Polarization 
Monostatic  RCS  at  6  GHz 
(CRAY  X-MP/216) 


numerical  results  reveals  a  pretty  good  correlation  between  the  two  when 
defining  the  troughs  and  the  nodes.  The  magnitude  of  the  RCS  face  on  (full 
side  of  disk)  is  in  good  agreement,  but  as  the  measurements  near  the  "thin" 
edge  the  similarity  falls  off.  Possibly  this  is  a  region  that  doesn't 
contain  enough  mesh  points  per  wavelength,  lambda,  and/or  a  good  enough 
aspect  ratio  to  get  a  good  converged  solution.  In  reference  [5],  user's  of 
EM-TRANAIR  at  the  Naval  Weapon  Center  at  China  Lake  had  problems  with  the 
numerical  solutions  on  closed  cylinders. 

"[They] .. .discovered  that  solutions  for  vertically  polarized  scattering  from 
conducting  cylinders  differed  significantly  from  experimental  measurements 
and  from  solutions  obtained  by  other  computational  techniques."  The 
conducting  circular  disk,  being  a  closed  cylinder  may  be  incurring  this 
inherent  problem.  More  field  testing  is  encouraged. 
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CONCLUSION 


> 

EM-TRANAIR  is  a  computer  program  for,the  solution  of  Maxwell's 
equations  in  three  dimensions.  This  primer  on  EM-TRANAIR  installation  at 
WRDC  and  guide  to  confronted  execution  problems  will  aid  in  understanding 
the  internals  of  EM-TRANAIR.  Successful  pathways  for  accurate  converged 
numerical  solutions  are  possible. 

The  architecture  of  EM-TRANAIR  was  briefly  presented  in  the  discussion 
of  various  installation  problems  encountered  when  moving  from  a  CFT  compiler 
to  a  CFT77  compiler.  Compilation  problems  and  solutions  deal!*  with  passing 
arguments,  dimensioning  arrays  and  boolean  algebra  manipulations. 

Execution  problems  with  the  input  processor,  V1INP  and  the  input 
solver,  V1S0L  were  presented.  Troubleshooting  methodologies  and  working 
solutions  .were  revealed.  The  actual  problems  needing  to  be  surmounted  for 
successful  program  execution  are  the  examples  contained  in  this  document. 

Different  input  configurations  were  used  to  establish  "rule-of-thumb" 
parameters  to  achieve  accurate  numerical  solutions.  Parameters  such  as 
incidence  frequency,  number  of  mesh  points  and  aspect  ratioswere  discussed. 

Finally,  the  area  to  prove  being  the  most  successful  in  the  post¬ 
processing  examination  was  the  Rayleigh  region. 

The  numerical  RCS  solutions  presented  are  proof  positive  that  the 
EM-TRANAIR  code  is  executable  and  providing  good  computational  RCS  results. 
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APPENDIX  A 


SUBROUTINE  LISTINGS  FOR  LIBRARIES,  INPUT  PROCESSOR  AND  SOLVER 

This  appendix  is  intended  to  guide  the  user  through  the  various 
pathways  as  the  result  of  CALL  SUBROUTINE  statements.  The  listings  were 
generated  with  UPDATE  and  the  CD  output  option  under  COS. 
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♦****APPENDIX  A 

-  V1INP 

SUBROUTINE  DECKS 

******** 

**************** 

***♦♦* 

AAINPUT 

INIT 

FLDFIL 

ERRMSG 

NETMOD 

OPDEF 

PIW4 

UNPIW4 

BOPDEF 

EDGPNT 

FBFND 

FESBLE 

FLDDEF 

LPCUT 

LTCUT 

PANDQC 

PANDEF 

PANPNT 

PFBLST 

PIFBCH 

TADJST 

TCOORD 

TIFBCH 

TRIPNT 

STRTAU 

ABTCHK 

CPETP 

EDGABT 

DETEQV 

PROEQV 

CHGGRD 

PRONEQ 

PEADET 

NETWRKS 

MSHIND 

EMARK 

FBXLST 

PADJST 

PANGHB 

EPTAU 

ABTIDN 

ABTDNC 

DNCHEK 

MPNORM 

D2LINE 

DOMVEC 

ABTAIP 

ABTDZO 

DZBCHK 

ABTECD 

ABTFSD 

ABTDAB 

ABTSYM 

ABTDAI 

ABTDUE 

ABTAIO 

ABTABO 

FSGCMP 

ADDIN2 

ABTEMP 

ABTEQC 

MPTEQC 

ABTMSG 

ABTEND 

CMPIED 

EDGIND 

EPOINT 

NRMESH 

ABTELS 

IDNGEO 

ABTPOS 

EDGMPI 

NREDGE 

INTNRM 

EDGTAU 

OUTPKV 

BMARK 

PERPNT 

FILOPN 

KEYINV 

OPGEN 

PEAFOR 

PROMIN 

LSOLVE 

MINCUT 

GRDGEN 

MSPNT1 

CAMBER 

CIRC 

ELLPT 

GADNET 

TRINP 

I  NBC 

INBC1 

INBC2 

MESHP 

QUADNT 

TRWAKE 

PEAINP 

MSHDNS 

INECHO 

INSUM 

IJKLAT 

LATIJK 

GRDDEF 

ASEGL 

ABTSEG 

APNTL 

ABTPTl 

ABTPT2 

CHCOLP 

BCONCL 

SPBSC 

BCOPT 

BLCAL 

BLCCAL 

CMNGRD 

CONTRL 

CPABT 

CPNOR 

DASPL 

XDASPL 

XDASIN 

EDGSPL 

GEOMC 

REFLOC 

SINCD 

SING 

SNGDEL 

SUBPQR 

SURPRO 

TCNTRL 

TSING 

ZCADJ 

ENRCHG 

GTALAM 

NRPTHP 

INSIDE 

NRPTED 

OPSTNC 

DEFOPG 

OPSET 

DEFDPT 

APTDEL 

DEFFPT 

TFILL 

ADJWK 

DREGD 

VOPCAL 

MAP  I NV 

EQUIVC 

LSCHOP 

EIFBCH 

CIFBCH 

TRICAL 

ECAL 

FCAL 

HCAL 

SOPCAL 

PDRVR 

GLOCAT 

VELCAL 

TBTOP 

PEQCL 

PINSID 

PNGHBR 

INNGBR 

EXNGBR 

REQC 

TREG 

PTINT 

ABTORD 

VMOMC 

EMNCUT 

CMPCTL 

OPPLOT 

DEFPPT 

SMOMC 

LINTRI 

INDCHK 

TRICPI 

FPNPLT 

SELUL 

GNRML 

FLDCHK 

TRINRM 

PAREA 

TESTP 

OPCHK 

MATCHK 

FAOPRT 

AGNRML 

CDRVR 

CNTD 

VCALF 

VCALA 

IBSET 

DEFCPT 

LINLS 

CDRVR1 

IBSET1 

DFCPT1 

VCALAA 

LUMPC 

LUMPL 

LUMPS 

LSQSFT 

BASEVL 

BASISQ 

ABTINT 

GPHSCN 

GPHPLK 

GPLUCK 

DERIVE 

SITPNT 

THNBCN 

TKCAL 

TREBRN 

TRERUT 

KID 

FACRDC 

FACTOP 

TDPEQV 

ODPF.QV 

ANGORD 

EDGRDC 

PLNEQV 

PROJD 

SOUT 

WRFAOP 

WRSTOP 

FREEOP 

BDYST 

WVCAL 

CAMLT2 

TZCAL 

CORIND 

FLDSCT 

FSDCOF 

HNGHBR 

MXMT 

OPCOMP 

DMATOP 

EMOUT 

OPADD 

LSTDEF 

CONDEF 

CJDEF 

SRCDEF 

RHODEF 

GRHDEF 

GAMDEF 

*♦*♦* APPENDIX  A 

-  GFLIB 

SUBROUTINE  DECKS 

************ 

******** 

********** 

INICOR 

SETCOR 

GETCOR 

FRECOR 

BNQUAD 

BVPGEN 

BVPG3D 

DFNG3D 

DFNIJK 

DFNSGV 

GETGVL 

GQXTND 

GRECUR 

OUTPKV 

GQRG12 

G3DABT 

GQRG3 

TIMING 

OMGCON 

GRNAXS 

COMGEN 

LSQSOL 
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LSQFBS 

RDG3D 

VIP 

ZERO 

DCBHTG 

OUTLIN 

OUTMAT 

OUTVEC 

PIW4 

UABEND 

UNPIW4 

AMCON 

SHELLX 

SPARSA 

SPARSE 

SPARSF 

SPARSO 

SPARSP 

SPARSS 

POLGEN 

ASYG3D 

APPENDIX  A 

-  EMLIB  : 

SUBROUTINE  DECKS*  **** • *  * 

*  *  *  *  *  *  * 

FFZER 

SPCSCAL 

SPDIV  SPGRAD 

DIV 

CCONT 

CCD ISP 

CFILLI 

FCOPY  CAXPYI 

CDOTCI 

CGTHRS 

WHENCGT 

VHEAD  CASSUM 

BUBBLE 

QUICK 

DIVIDE 

CCDSPY 

APPENDIX  A 

-  GPLIB 

SUBROUTINE  DECKS 

********** 

******************** 

CADD 

CFILL 

CGATHR 

CSCATT 

CZERO 

RSIGN 

ICONTR 

CCONTUR 

CONTUR 

CSTPRT 

CVECWR 

DSP 

DSPPLT 

CONJUG 

RMATWR 

IMATWR 

OVECWR 

CMATWR 

PLTCAL 

PPLOT 

FSHELL 

KEYSRT 

7SHELL 

RELIST 

SETPT2 

SETPT1 

CPLOT 

ICOPY 

I  ZERO 

LSPARI 

PARINT 

RFILL 

RZERO 

LISTN 

IADD 

NRMARG 

CRXIPY 

CRXPYI 

CXITY 

OMATWR 

RTENWR 

CTENWR 

ITENWR 

OTENWR 

CROSS 

DISTNC 

MAG 

MNMOD 

PIDENT 

SHLSRT 

UKYSRT 

UVECT 

VADD-CR 

VMUL-CR 

IFILL 

CCXIPY 

CCXPYI 

MODVEC 

PlwJrtJ 

PRCONT 

ANSPEK 

CGMRES 

GMRES 

NGMRES 

aiTRXTP 

UNIONL 

INTRSL 

SETPTR 

se: .HR 

CMXM 

CMXMA 

RRXIPY 

RRXPYI 

RXITY 

BSEQR 

BSEQW 

REVERS 

GWRANA 

MSNDX 

CQLSF 

DCBHTX 

LSQSFX 

DEFVEC 

DCBHT 

SGINV 

LSQSF 

I'SIjv  J'P 

LSQSG 

INCMPR 

SCMPKT 

SORTAK 

TCOF 

TRANS 

SRCHOL 

CPLOT1 

MTVCWR 

NORCAL 

IBSRCH 

UVECTM 

BINSCH 

INTP1 

INTP2 

QCOF 

NFILL 

BISEQF 

BISEQW 

UMAX 

IIMIN 

PCCONT 

CHCONT 

ROUND 

RREWND 

BISEQR 

HSMCON 

CMUL 

HSFFTS 

HSFFT1 

HSFFT2 

HSFFT3 

HSFFT4 

HSFFT5 

HSFFT6 

HSFFT7 

HGERR 

MSCAL 

MSUM 

MCOPY 

M2COPY 

MAXPY 

M2AXPY 

M2SCAL 

MSCOPY 

M2SCPY 

M2SBTR 

MSBUTR 

MBUTR 

M2SBFT 

MPSUM 

M2PSUM 

MVADD 

MSWAP 

H1FFTR 

H2FFTR 

MARITH 

H1SIN 

H1SSIN 

H1COS 

PERM2P 

RCMUL 

PERM2 

PERMS 

PERM2S 

M2ZERO 

MVADSB 

MZERO 

FACTOR 

MPERM 

M2PERM 

HFFTS 

BFFTP 

H2SIN 

H2C0S 

H2SSIN 

MVCMUL 

MADSB 

VMCOPY 

MOUT 

MI  OUT 

PGOUT 

HlSCOS 

H2SCOS 

SSCOPY 

ROPEN 

RREAD 

RWRIT 

RCLOS 

MSHIFT 

MSSUM 

GETLIN 

RTSRCH 

TSRCH 

RVECWR 

IVECWR 

PYCONT 

ISMX64 

SGBSL 

MASUM 

MFILL 

MVSADD 

SGBL 

SGBLT 

SGBU 

SGBUT 

INDFND 

PYCSRT 

CPMAT 

IMSHFT 

IMFILL 

VN0RM2 

REINT 

SGBFA 

SGBCO 

SETRA 

GETRA 

IOFRA 

GCOL 

FRERA 

SETCIO 

FNDREC 

REDRA 

GACRA 

ISCATC 

ISCATZ 

CNTIP 

CNTRP 

MAPEQV 

LSRUT 

TREBRN 

TRERUT 

KID 
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CCALN 

COMPIP 

MSROTM 

PANUNI 

PROJ 

PANMOM 

ROT2 

SCALE 

UNIPAN 

POILO 

H1SYM 

H2SYM 

GBPRO 

GBPRY 

GRADML 

GYADML 

GREEN 

REGGRN 

BDYCVO 

BADD 

EXGRN 

PERMZ 

POIFO 

POIFY 

CMPSCL 

GCPCAL 

SURFIT 

XXADJ 

LCHVAR 

EDGLS 

TRAN 

ZMPROJ 

POILY 

P0IL1 

BGMUL 

BGMUL 1 

GBPRYZ 

GDB 

GDBO 

GRNASM 

GREEN 1 

SMGRN 

BCOP 

BADDO 

POISOL 

EXP OIO 

POIF1 

THRMIS 

GRDIND 

ISCAL 

KRMOVE 

LPROJ 

ROTl 

ZMERGE 

TRI2D 

TRI2DS 

GRNDWN 

GRPLAN 

GYPLAN 

GCOMP 

BDYCV1 

BXSIZE 

BDYCVY 

EXPOIY 

EXPOIl 

*****APPENDIX  A  -  EXSOL  SUBROUTINE 


1/UVI\U  ' 


H1SIN 

H2SIN 

H1FFTR 

H2FFTR 

PERM 

PERM1P 

PERM2P 

PERMS 

MADD 

SWTSUM 

TRI3H 

TRI2H 

TRI2R 

PERMC 

H1C0S 

H2C0S 

CONBOX 

CONVLU 

CONFFT 

GBPRHE 

MVADDN 

HELMOP 

HELMO 

EXHELM 

BDYCVH 

HCFFTS 

HCFFT1 

HCFFT2 

HCFFT3 

HCFFT4 

HCFFT5 

HCFFT6 

HCFFT7 

MSHUFL 

M2SHFL 

ISIN11 

HFIRST 

HLAST 

GTRANS 

MABUTR 

» * » *  * APPENDIX  a  -  V1SOL  SUBROUTINE  DECKS **♦**». *********************** 


AA3DS 

ACAL 

ALTFSP 

AQCAL 

BLDNDX 

BLDIEL 

CGMRES 

DDIVG 

DERIVE 

DGRAD 

DJCOLL 

DNSADD 

DNSTFL 

FLTDNS 

FOPEN 

GREENF 

HSOLVE 

INCFAD 

INICOM 

INIT 

JINIT 

JSTORE 

JWRITE 

OQCAL 

OUTPUT 

OCVWR 

READIN 

RDFAOP 

SETUP 

SOLVER 

TCAL 

WVCAL 

CLUSS 

CLUSM 

LCAL 

LICAL 

CLUSO 

CLUSB 

BCCAL 

OPGET 

RADSUB 

MDBIAS 

ORDER 

PROXIM 

RADBC 

CORIND 

REDSET 

ADJAC 

OPCOUP 

INFOLV 

CMPTLW 

FNDLV 

BLDLW 

EQUIVC 

MODVEC 

FFRCS 

BISRCS 

GAP ARM 

XFMMTX 

MONRCS 

CJGDJ 

MONOUT 

EQNRDC 

FNDIND 

INDFND 

INDMRG 

INDSTD 

RADELM 

RADFND 

RADIND 

RDCEQN 

REORDR 

RESTAK 

RSINDX 

PSIELM 

LINORD 

MOMEQ 

PLAORD 

PLNORD 

NESTED 

SELPLN 

PSHSTK 

REFORM 

CSINIT 

CSINPT 

CSFLTR 

CSFLSH 

CSPREP 

CSSORT 

CSSRTB 

CSRADX 

CSMRGO 

CSMRGI 

CSADDO 

CSADDI 

CSCNRC 

CSTRNF 

CSLU 

CSDEC 

DMPRB 

FCHRB 

INIRB 

STORB 

SLSTNE 

LLSTEQ 

SGTHR 

CGTHR 

CGTHRZ 

CSCTR 

ROWSCL 

INCFDI 

IEQUAL 
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APPENDIX  B 


ERROR  MESSAGES  GENERATED  BY  THE  CFT77  COMPILER  ON  V1INP 

This  appendix  contains  a  listing  of  the  error  messages  that  were 
returned  when  the  original  input  processor  source  deck  was  compiled  using 
the  FORTRAN  CFT77  compiler.  The  code  was  compiled  using  a  FORTRAN  compiler 
CFT  when  it  was  developed  by  Boeing.  Error  messages  will  lead  to  fatal 
crashes  when  decks  containing  them  are  tried  to  execute.  Warning  messages 
may  lead  to  fatal  crashes  and  cautions  are  usually  just  noteworthy  comments. 

Use  this  listing  or  others  that  can  be  generated  for  the  library  decks 
and  the  solver  code  setting  ON“H  in  the  CFT77  job  control  statement  under 
COS.  These  can  be  compared  to  the  successful  executable  decks  to  trace  down 
any  changes  the  author  made  to  the  EM-TRANAIR  code. 
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APPENDIX  C 


ORIGINAL  AND  REVISED  INITIALIZATION  FOR  VLINP 

The  two  decks  listed  in  this  appendix  were  used  to  initialize  variables 
for  the  input  processor  V1INP  source  program  AAINPUT.  The  first  listing  is 
from  the  original  code  where  a  call  statement  in  AAINPUT  is  used  to  call  up 
the  BLOCK  DATA  INIT.  The  second  listing  is  the  result  of  solving  an 
execution  error  by  imbedding  the  initialization  in  the  source  program 
AAINPUT. 
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♦  * » *  *  APPENDIX  C  -  INITIALIZATION  OF  PARAMETERS  FROM  ORIGINAL  INPUT 
PROCESSOR  DECK  V1INP* ♦♦♦♦♦♦♦♦♦**♦*******»****»**♦***♦************* 

♦DECK  INIT 

BLOCK  DATA  INIT 
♦CA  PARAM 
♦CA  /CHRCTR/ 

♦CALL  PRNT 
♦CALL  ABTPRT 
♦CALL  ABTNEW 
♦CA  /CRSDAT/ 

♦CA  /FILDAT/ 

♦CA  /GRIDQ/ 

♦CA  /INCFEM/ 

♦CA  /ITERQ/ 

♦CALL  /ARV1S/ 

♦CA  /NUMRCN/ 

♦CA  /PHYSCN/ 

♦CA  /LSCALE/ 

♦CALL  ACASE 
♦CALL  COMPRS 
♦CALL  SYMM 
♦CALL  /OFBOD/ 

♦CALL  KSTMLN 
♦CALL  FMCOF 
C 

C  INITIALIZE  FUNDAMENTAL  COMMON  BLOCK  VARIABLES 
C 

C  /CHRCTR/ 

DATA  REVDAT/ ’7/31/87 ’/ 

DATA  TITLE  /’  DEFAULT  TITLE 

1  7,  NAME  /’  DEFAULT  NAME 

2  7, 

3  IPW  /’PW  7,  IPOL  /  ’EFLD  ’/ 

C  /PRNT/ 

DATA  IGEOMP  /O/,  ISINGP  / 0/ ,  ICONTP  /O/,  IBCONP  /O/ 

DATA  IEDGEP  fO/ ,  ISINGS  /O/,  IPRAIC  /O/,  IPARTP  /O/ 

DATA  IPARTS  /O/ ,  NEXDGN  /Of,  IOUTPR  /O/,  IFMCPR  /O/ 

DATA  ICOSTP  /O/ ,  IEXTRP  /O/,  ISPMAP  /O/,  ICPMAP  /O/ 

DATA  IBCMAP  /O/ 

C  /ABTPRT/ 

DATA  I GEOIN  /!/,  IGEOUT  /O/,  NWXREF  /O/,  NWPROP  /O/,  IABUTD  /O/ 
DATA  IABSUM  /l/,  IFABPR  /l/ 

DATA  IGEOM  /O/,  IBCND  /O/,  IOPER  /0/t  IPANL  /Of 
DATA  NWDK  /O/ 

C  /ABTNEW/ 

DATA  EPSGEO  /O./ 

C  /ACASE/ 

DATA  ALPHA  /4»0./,  BETA  /4^0./,  FSVM  /4*1 ./ 

DATA  I ACASE  /!/,  NACASE  /l/ 

C  /COMPRS/ 

DATA  AMACH  /0.0001/,  ALPC  /<)./,  BETC  /07 
C  /SYMM/ 

DATA  NSYMM  /O/,  MISYM  /O/,  MJSYM  /O/ 

C  /OFBOD/ 

DATA  NOF  /Of 
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C  /KSTMLN/ 

DATA  NUMPTS  /Of 
C  /FMCOF/ 

DATA  SREF  /I ./,  BEEF  /l-/,  CREF  /l./(  DREF  /I ./ 

DATA  XREF  /O./,  YREF  /O./,  ZREF  /O. / 

C  /FILDAT/ 

DATA  NIN  /5/ ,  NOUT  /6/ ,  NPRSL  /7/,  NPRPO  /8/ 

DATA  NABFL  /1 1/ ,  NFAKI  /12/,  NFBLP  /13/,  NFBKP  /14/, 

NFVD / 1 5/ , NFVC/ 1 6/ , NFBC/ 1 7/ , NFGRE / 1 8/ , NFGRD / 1 9  ' , 

NFDWF/20/ , NFDWA/21/ , NFDWD/22/ , NFPTL/23/ , NFMTL/24/ , 
NFSCR/25/ 

C  /CRSDAT/ 

DATA  NABSW  /l/,  NANGB(l)  /l/,  ANGB(1,1)  /O./,  ANGB(2,1)  /O .  / , 
2  ANGB (3,1)  /O . / 

C  /GRIDQ/ 

DATA  NX  /NGRDl/ ,  NY  /NGRD2/ ,  NZ  /NGRD3/ 

DATA  XI  /O. /,  XF  10./,  YI  /O./,  YF  /O./,  ZI  /O . / ,  ZF  /O . / 

C  /INCFEM/ 

DATA  FREQ/3 . E8/ , EHMAG / ( 1 . , 0 . ) / , 

2  NASW  /If,  NANG ( 1 )  /l/,  ANG(l.l)  /O . / ,  ANG(2,1)  /O./, 

3  ANG(3, 1)  /O. / 

C  /ITERQ/ 

DATA  TOL  /l.E-10/,  NITER  /SO/,  NSRCH  /50/,  MSGLVL  /O/ 

DATA  RLOSS  /O./,  DROPT  /O . / ,  NRELXD  /l/ 

DATA  IDMRCH  /3,5*0/,  IDSWP  /2,5*0/ 

C  /ARVIS/ 

DATA  NVIS  /l/,  RMCUT  / . 95/ 

C  /LSCALE/ 

DATA  RLEN  /I .  / 

C  /NUMRCN/ 

DATA  PI  /3. 141592653358979/,  CONE  / (0. , 1 . ) / 

C  /PHYSCN/ 

DATA  GAMMA  /I .4/ 

DATA  VLIGHT/2 . 997925E8/ 

END 
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♦♦♦♦♦APPENDIX  C  -  INITIALIZATION  IMBEDDED  IN  INPUT  PROCESSOR  V1INP 
SOURCE  PROGRAM  AAINPUT* ******************************************* 

.DECK  AAINPUT 

PROGRAM  AAINPUT 
C 

C  DRIVER  FOR  THE  TRANAIR  INPUT  PROCESSOR 
C 

*CA  /CASE/ 

*CA  /FILDAT/ 

.CALL  ABTPRT 
.CALL  /CHRCTR/ 

.CALL  /GRIDQ/ 

*CA  PARAM 
.CALL  PRNT 
.CALL  ABTNEW 
*CA  /CRSDAT/ 

*CA  /INCFEM/ 

*CA  /ITERQ/ 

.CALL  /ARVIS/ 

*CA  /NUMRCN/ 

*CA  /PHYSCN/ 

*CA  /LSCALE/ 

.CALL  ACASE 
.CALL  COMPRS 
.CALL  SYMM 
.CALL  /OFBOD/ 

.CALL  KSTMLN 
.CALL  FMCOF 

REAL  CARD (10) 

C 

C  INITIALIZE  FUNDAMENTAL  COMMON  BLOCK  VARIABLES 

C 

C  /CHRCTR/ 

DATA  REVDAT/ ’7/31/87 7 
DATA  TITLE  /’  DEFAULT  TITLE 

1  7,  NAME  /’  DEFAULT  NAME 

2  7, 

3  IPW  /  ’PW  7,  IPOL  /’EFLD7 

C  /PRNT/ 

DATA  IGEOMP  /O/ ,  ISINGP  /O/ ,  ICONTP  /O/,  IBCONP  /O / 

DATA  IEDGEP  /O/ ,  ISINGS  /O/ ,  IPRAIC  /O/,  IPARTP  /O/ 

DATA  IPARTS  /Of,  NEXDGN  /O/,  IOUTPR  /O/,  IFMCPR  /O/ 

DATA  ICOSTP  /O/,  IEXTRP  /O/ ,  ISPMAP  /O/,  ICPMAP  /O/ 

DATA  IBCMAP  /O/ 

C  /ABTPRT/ 

DATA  I GEOIN  /!/ ,  IGEOUT  /Of,  NWXREF  /O/,  NWPROP  /O/,  IABUTD  /O/ 
DATA  IABSUM  /l/,  IFABPR  /l/ 

DATA  IGEOM  /O/,  IBCND  /O/,  IOPER  /O/,  IPANL  /Of 
DATA  NWDK  /O/ 

C  /ABTNEW/ 

DATA  EPSGEO  /O./ 

C  /ACASE/ 

DATA  ALPHA  /4*07,  BETA  /A*0./,  FSVM  /4*17 
DATA  IACASE  /l/,  NACASE  /l/ 

C  /COMPRS/ 
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DATA  AMACH  /0.0001/,  ALPC  /O./,  BETC  /O .  / 

C  /SYMM/ 

DATA  NSYMM  /O/,  MISYM  /O/,  MJSYM  /O/ 

C  /OFBOD/ 

DATA  NOF  /O/ 

C  /KSTMLN/ 

DATA  NUMPTS  JO/ 

C  /FMCOF/ 

DATA  SREF  Jl.J,  BREF  /!./,  CREF  /!./,  DREF  Jl.J 
DATA  XREF  JO./,  YREF  JO./,  ZREF  JO./ 

C  /FILDAT/ 

DATA  NIN  /5/ ,  NOUT  /6/,  NPRSL  J7J,  NPRPO  /8 / 

DATA  NABFL  /ll/,  NFAKI  /12/,  NFBLP  /13/,  NFBKP  /14/, 

NFVD/15/ .NFVC/16/ , NFBC/17/ ,NFGRE/18/ ,NFGRD/19/ , 

NFDWF/20 / , NFDWA/21 / , NFDWD /22 / , NFPTL/23/ , NFMTL/24/ , 
NFSCR/25/ 

C  /CRSDAT/ 

DATA  NABSW  /I/,  NANGB(l)  /I/,  ANGB(1,1)  JO./,  ANGB(2,1)  JO./, 
2  ANGB(3, 1)  JO./ 

C  /GRIDQ/ 

DATA  NX  /NGRD1 / ,  NY  /NGRD2/,  NZ  /NGRD3/ 

DATA  XI  JO./,  XF  JO./,  YI  JO./,  YF  JO./,  ZI  JO./,  ZF  JO./ 

C  /INCFEM/ 

DATA  FREQ/3. E8/,EHMAG/(1. ,0.)/, 

2  NASW  /!/ ,  NANG(l)  J\J ,  ANG(1,1)  JO./,  ANG(2,1)  JO./, 

3  ANG(3 , 1)  JO./ 

C  /ITERQ/ 

DATA  TOL  /l.E-10/,  NITER  /50/ ,  NSRCH  /50/,  MSGLVL  JO/ 

DATA  RLOSS  JO./,  DROPT  JO./,  NRELXD  /!/ 

DATA  IDMRCH  / 3,5*0/,  IDSWP  /2,5*0/ 

C  /ARVIS/ 

DATA  NVIS  /If,  RMCUT  / .95/ 

C  /LSCALE/ 

DATA  RLEN  /I . / 

C  /NUMRCN/ 

DATA  PI  /3. 141592653358979/,  CONE  / (0. , 1 . ) / 

C  /PHYSCN/ 

DATA  GAMMA  /I .4/ 

DATA  VLIGHT/2 . 997925E8/ 

C 

WRITE (NOUT, 6005) 

C 

CALL  DERIVE 
C 

CALL  VHEAD(NOUT, ’V1INP: ’.REVDAT) 

WRITE (NOUT, 6030) 

C 

C 

CALL  CSTPRT(’ INITIAL’) 

C 

IERR=0 

CALL  TRINP(IERR) 

C 

CALL  CSTPRT(’TRINP’) 

IF(IERR.EQ.3)  THEN 
WRITE (NOUT, 6060) 
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CALL  ABORT (’END-OF-FILE  ON  INPUT’) 

ENDIF 

C 

IF (IERR.NE.O)  THEN 

WRITE (NOUT, 6050)  I ERR 

CALL  ABORT 

ENDIF 

CALL  OPDEF 

CALL  CSTPRT( ’OPDEF’) 

DO  130  1=1,4 

IF (IDCHK(I) .EQ. -1)  GO  TO  200 

130  CONTINUE 
C 

CALL  SOUT 

CALL  CSTPRT ( ’SOUT ’ ) 

C 

200  CONTINUE 
C 

READ(NIN,END=210)  card 
GOTO  200 
C 

210  CONTINUE 
C 

REWIND  NPRSL 
REWIND  NPRPO 
C 
C 

WRITE (NOUT, 6040) 

C 

C 

6005  FORMAT(///10X, ’ENTERING  INITIALIZATION  PHASE’) 

6030  FORMAT (///l OX, ’  STARTING  INPUT  PROCESSOR  PROGRAM’) 
6040  FORMAT(//10X, ’COMPLETING  INPUT  PROCESSOR  PROGRAM’) 
6050  FORMAT <//l OX, ’PROGRAM  ABORT  DUE  TO  ERROR’, 14) 

6060  FORMAT(//10X, ’END-OF-FILE  ON  INPUT  FILE’) 

STOP 

END 
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APPENDIX  D 


ORIGINAL  AND  REVISED  CONDUCTING  SPHERE  NETWORK  INPUT  DATA 

This  appendix  contains  the  original  and  revised  listings  for  the 
conducting  sphere  configuration.  The  original  network  configuration  used  by 
Boeing  contains  twice  as  many  points  as  the  executable  network. 
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♦♦.♦♦APPENDIX  D  -  ORIGINAL  CONDUCTING  SPHERE  SURFACE  GRID  NETWORK  INPUT 
DATA  *****♦♦♦♦  +  ♦♦♦♦***♦♦♦*♦*♦**♦♦  +  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦**♦♦♦♦♦♦♦♦♦♦♦♦♦ 


STITLE 
SPHERE 
CASE  #1 
SSYM 
0. 

SEAT 


NO  PLANES  OF  SYMMETRY 


0. 


0.0 

-1 . 

0.0 

0.0 

0.0 

SDAT 

0. 

♦  1  . 

0. 

0. 

0. 

0.0 

SBOX 

- 1 . 0000 

1.0000 

17. 

-1 .0000 

1 .0000 

17. 

- 1 . 0000 

1.0000 

17. 

SNET 

SCIR 

A 

4  . 

1 . 

1. 

0.0 

0. 

RIGHT 

41  . 

.8 

0. 

.79938323 

.03140785 

.79753387 

.06276728 

.79445477 

.09402992 

.79015067 

.12514757 

.78462822 

.15607226 

.77789594 

. 18675629 

.76996419 

.21715236 

.76084521 

.24721360 

. 75055307 

.27689365 

.73910363 

.30614675 

.72651454 

.33492779 

.71280522 

.36319240 

.69799681 

.  39089699 

.68211213 

.41799885 

.66517569 

.44445619 

.64721360 

.47022820 

.62825355 

.49527516 

.60832477 

.51955844 

.58745801 

.  54304060 

.56568543 

.56568543 

. 54304060 

.58745801 

.51955844 

. 60832477 

.49527516 

.62825355 

. 47022820 

.64721360 

.44445619 

.66517569 

.41799885 

.68211213 

.39089699 

.69799681 

.36319240 

.71280522 

.33492779 

.7265145 4 

.30614675 

.73910363 

.27689365 

. 75055307 

.24721360 

. 76084521 

.21715236 

.76996419 

.  18675629 

.77789594 

.15607226 

.78462822 

. 12514757 

.79015067 

.09402992 

.79445477 

.06276728 

. 79753387 

.03140785 

.79938323 

0.0 

.8 

81  . 

90. 

92.25 

94.5 

96.75 

99.0 

101.25 

103.5 

105.75 

108. 

110.25 

112.5 

114.75 

117. 

119.25 

121.5 

123.75 

126. 

128.25 

130.5 

132.75 

135.0 

137.25 

139.5 

141.75 

144. 

146.25 

148.5 

150.75 

153. 

155.25 

157.5 

159.75 

162.0 

164.25 

166.5 

168.75 

171. 

173.25 

175.5 

177.75 

180. 

182.25 

184.5 

186.75 

189. 

191.25 

193.5 

195.75 

198. 

200.25 

202.5 

204 . 75 

207.0 

209.25 

211.5 

213.75 

216.0 

218.25 

220.5 

222.75 

225. 

227.25 

229.5 

231.75 

234.0 

236 . 25 

238.5 

240.75 

243.0 

245.25 

247.5 

249.75 

252.0 

254.25 

256.5 

258.75 

261 .0 

263.25 

265.5 

267.75 

270.0 

0. 

41 

0. 


LEFT  FRONT 


.8 


-.03140785.79938323 


- . 09402992 . 79445477  - . 12514757 . 79015067 


-.06276728.79753387 

-.15607226.78462822 
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-.18675629.77789594  -.21715236.76996419  -.24721360.76084521 
-.27689365.75055307  -.30614675.73910363  -.33492779.72651454 
-.36319240.71280522  -.39089699.69799681  -.41799885.68211213 
-.44445619.66517569  -.47022820.64721360  -.49527516.62825355 
-.51955844.60832477  -.54304060.58745801  -.56568543.56568543 
-.58745801.54304060  -.60832477.51955844  -.62825355.49527516 
-.64721360.47022820  -.66517569.44445619  -.68211213.41799885 
-.69799681.39089699  -.71280522.36319240  -.72651454.33492779 
-.73910363.30614675  -.75055307.27689365  -.76084521.24721360 
-.76996419.21715236  -.77789594.18675629  -.78462822.15607226 
-.79015067.12514757  -.79445477.09402992  -.79753387.06276728 


-.79938323.03140785 

O  1 

-.8 

0.0 

Ol  . 

90. 

92.25 

94.5 

96.75 

99.0 

101.25 

103.5 

105.75 

108. 

110.25 

112.5 

114.75 

117. 

119.25 

121.5 

123.75 

126. 

128.25 

130.5 

132.75 

135.0 

137.25 

139.5 

141.75 

144. 

146.25 

148.5 

150.75 

153. 

155.25 

157.5 

159.75 

162.0 

164.25 

166.5 

168.75 

171. 

175.25 

175.5 

177.75 

180. 

182.25 

184.5 

186.75 

189. 

191.25 

193.5 

195.75 

198. 

200.25 

202.5 

204 . 75 

207.0 

209.25 

211.5 

213.75 

216.0 

218.25 

220.5 

222.75 

225. 

227.25 

229.5 

231.75 

234.0 

236.25 

238.5 

240.75 

243.0 

245.25 

247.5 

249.75 

252. 

265.5 

0. 

41. 

254.25 

267.75 

256.5 

270.0 

258 . 75 

261.0 

263.25 

RIGHT 

.8 

0. 

.79938323 

.03140785 

.79753387 

.06276728 

.79445477 

.09402992 

.79015067 

.12514757 

.78462822 

.15607226 

.77789594 

.18675629 

.76996419 

.21715236 

.76084521 

.24721360 

. 75055307 

.27689365 

.73910363 

.30614675 

.72651454 

.33492779 

.71280522 

.36319240 

.69799681 

.39089699 

.68211213 

.41799885 

.66517569 

.44445619 

.64721360 

.47022820 

.62825355 

.49527516 

. 60832477 

.51955844 

.58745801 

. 54304060 

.56568543 

56568543 

. 54304060 

.58745801 

.51955844 

.  60832477 

.49527516 

. 62825355 

. 47022820 

.64721360 

.44445619 

.66517569 

.41799885 

.68211213 

. 39089699 

.69799681 

.36319240 

.71280522 

.33492779 

.72651454 

.30614675 

.73910363 

.27689365 

.75055307 

.24721360 

.76084521 

.21715236 

.76996419 

.18675629 

.77789594 

.15607226 

.78462822 

.12514757 

.03140785 

81. 

270. 

.79015067 

.79938323 

.09402992 

0.0 

.79445477 

.8 

.06276728 

.79753387 

272.25 

274.5 

276.75 

279. 

281.25 

283.5 

285.75 

288. 

290.25 

292.5 

294.75 

297. 

299 . 25 

301.5 

303 . 75 

306.0 

308.25 

310.5 

312.75 

315. 

317.25 

319.5 

321.75 

324. 

326 . 25 

328.5 

330.75 

333. 

335.25 

337.5 

339.75 

342. 

344.25 

346.5 

348.75 

351 . 

353.25 

355.5 

357.75 

000.0 

002.25 

004.5 

006.75 

009.0 

011.25 

013.5 

015.75 

018.0 

020.25 

022.5 

024.75 

027.0 

029 . 25 

031.5 

033.75 

036.0 

038.25 

040.5 

042.75 

045.0 

047.25 

049.5 

051.75 

054.0 

056.25 

058.5 

060.75 

063.0 

065.25 

067.5 

069.75 
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078.75 


83.25 


072.0  074.25  076.5 

85.5  87.75  90. 

0. 


81 . 


LEFT  BACK 


41 . 

0.  .8  -.03140785.79938323  -.06276728.79753387 

-.09402992.79445477  -.12514757.79015067  -.15607226.78462822 
-.18675629.77789594  -.21715236.76996419  -.24721360.76084521 
-.27689365.75055307  -.30614675.73910363  -.33492779.72651454 
-.36319240.71280522  -.39089699.69799681  -.41799885.68211213 
-.44445619.66517569  -.47022820.64721360  -.49527516.62825355 
-.51955844.60832477  -.54304060.58745801  -.56568543.56568543 
-.58745801.54304060  -.60832477.51955844  -.62825355.49527516 
-.64721360.47022820  -.66517569.44445619  -.68211213.41799885 
-.69799681.39089699  -.71280522.36319240  -.72651454.33492779 
-.73910363.30614675  -.75055307.27689365  -.76084521.24721360 
-.76996419.21715236  -.77789594.18675629  -.78462822.15607226 
-.79015067.12514757  -.79445477.09402992  -.79753387.06276728 
-.79938323.03140785  -.8  0.0 

81. 


270. 

272.25 

274.5 

276.75 

279. 

281.25 

283.5 

285.75 

288. 

290.25 

292.5 

294.75 

297. 

299.25 

301.5 

303.75 

306.0 

308 . 25 

310.5 

312.75 

315. 

317.25 

319.5 

321.75 

324. 

326.25 

328.5 

330.75 

333. 

335.25 

337.5 

339.75 

342. 

344.25 

346.5 

348.75 

351. 

353.25 

355.5 

357.75 

000.0 

002.25 

004.5 

006.75 

009.0 

011.25 

013.5 

015.75 

018.0 

020.25 

022.5 

024 . 75 

027.0 

029 . 25 

031.5 

033.75 

036.0 

038.25 

040.5 

042 . 75 

045.0 

047.25 

049.5 

051.75 

054.0 

056.25 

058.5 

060.75 

063.0 

065 . 25 

067.5 

069.75 

072.0 

074.25 

076.5 

078.75 

81  .0 

83.25 

85.5 

87.75 

90.0 

SITE 

50. 

5. 

0. 

0. 

50. 

1. 

0. 

0. 

0. 

SINC 

HFLD 

PW 

1.7893 

1. 

8. 

1. 

1 . 

0. 

0. 

0. 

0. 

0. 

0. 

SRCS 

2. 

91. 

0. 

0. 

0. 

2. 

0. 

0. 

91  . 

0. 

-90. 

90. 

2. 

0. 

0. 

STHE 

4. 

1. 

1 . 

2. 

0 

2. 

1. 

2. 

0 

3. 

1 . 

2. 

0 

4. 

1. 

2. 

0 
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.♦♦.♦APPENDIX  D  -  EXECUTABLE  CONDUCTING  SPHERE  NETWORK  WITH  FEWER 
SURFACE  GRID  POINTS* ********************************************* 

STITLE 

SPHERE  -  NO  PLANES  OF  SYMMETRY 

CASE  #1 

SSYM 

0.  0. 

SEAT 


0.0 

0.0 

0.0 

0.0 

0.0 

SDAT 

0. 

0.0 

0. 

0. 

0. 

0.0 

SBOX 

-1.0000 

1 .0000 

17. 

-1.0000 

1 .0000 

17. 

-1 .0000 
SNET 

1 .0000 

17. 

4. 

SCIR 


4. 

1 . 

0. 

41  . 

.8 

1 . 

0.0 

RIGHT  FRONT 

0. 

.79938323 

.03140785 

.79753387 

.06276728 

.79445477 

. 09402992 

.79015067 

.12514757 

.78462822 

.15607226 

.77789594 

.18675629 

.76996419 

.21715236 

.76084521 

.24721360 

. 75055307 

. 27689365 

. 73910 363 

.30614675 

.72651454 

. 33492779 

.71280522 

.36319240 

.69799681 

. 39089699 

.68211213 

.41799885 

.66517569 

.44445619 

.64721360 

. 47022820 

.62825355 

.49527516 

.60832477 

.51955844 

.58745801 

. 54304060 

.56568543 

.56568543 

.54304060 

.58745801 

.51955844 

.60832477 

.49527516 

.62825355 

.47022820 

.64721360 

.44445619 

.66517569 

.41799885 

.68211213 

. 39089699 

.69799681 

.36319240 

.71280522 

.33492779 

.72651454 

.30614675 

.73910363 

.2/689365 

.75055307 

.24721360 

.76084521 

.21715236 

.76996419 

.  1867o629 

.77789594 

.  15607226 

.78462822 

.12514757 

.03140785 

41. 

90. 

.79015067 
. 79938323 

.09402992 

0.0 

.79445477 

.8 

.06276728 

.79753387 

94.5 

99.0 

103.5 

108. 

112.5 

117. 

121.5 

126. 

130.5 

135.0 

139.5 

144. 

148.5 

153. 

157.5 

162.0 

166.5 

171 . 

175.5 

180. 

184.5 

189. 

193.5 

198. 

202.5 

207.0 

211.5 

216.0 

220.5 

225. 

229.5 

234.0 

238.5 

243.0 

247.5 

252  0 

0. 

41  . 

0. 

256.5 

261.0 

265.5 

270.0 

LEFT  FRONT 

.8 

- .0314078S 

..79938323 

-.06276728 

; . 79753387 

-.09402992.79445477  -.12514757.79015067  -.15607226.78462822 
-.18675629.77789594  -.21715236.76996419  -.24721360.76084521 
-.27689365.75055307  -.30614675.73910363  -.33492779.72651454 
-.36319240.71280522  -.39089699.69799681  -.41799885.68211213 
-.44445619.66517569  -.47022820.64721360  -.49527516.62825355 
-.51955844.60832477  -.54304060.58745801  -.56568543.56568543 
-.58745801.54304060  -.60832477.51955844  -.62825355.49527516 
-.64721360.47022820  -.66517569.44445619  -.68211213.41799885 
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-.69799681.39089699  -.71280522.36319240  -.72651454.33492779 
-.73910363.30614675  -.75055307.27689365  -.76084521.24721360 
-.76996419.21715236  -.77789594.18675629  -.78462822.15607226 
-.79015067.12514757  -.79445477.09402992  -.79753387.06276728 
-.79938323.03140785  -.8  0.0 

41 . 


90. 

94.5 

99.0 

103.5 

108. 

112.5 

117. 

121.5 

126. 

130.5 

135.0 

139.5 

144. 

148.5 

153. 

157.5 

162.0 

166.5 

171  . 

175.5 

180. 

184.5 

189. 

193.5 

198. 

202.5 

207.0 

211.5 

216.0 

220.5 

225. 

229.5 

234.0 

238.5 

243.0 

247.5 

252.0 

0. 

41. 

.8 

256.5 

261.0 

265.5 

270.0 

RIGHT  BACK 

0. 

.79938323 

.03140785 

.79753387 

.06276728 

. 79445477 

.09402992 

.79015067 

.12514757 

.78462822 

. 15607226 

. 77789594 

. 18675629 

.76996419 

.21715236 

.76084521 

.24721360 

.75055307 

.27689365 

.73910363 

.30614675 

.72651454 

.33492779 

.71280522 

.36319240 

.69799681 

.39089699 

.68211213 

.41799885 

.66517569 

.44445619 

.64721360 

. 47022820 

. 62825355 

.49527516 

. 60832477 

.51955844 

.58745801 

. 54304060 

.56568543 

.56568543 

. 54304060 

.58745801 

.51955844 

.60832477 

.49527516 

.62825355 

.47022820 

.64721360 

.44445619 

.66517569 

.41799885 

.68211213 

. 39089699 

.69799681 

.36319240 

.71280522 

.33492779 

.72651454 

.30614675 

.73910363 

.27689365 

. 7 5055307 

.24721360 

.76084521 

.21715236 

.76996419 

.18675629 

.77789594 

. 15607226 

.78462822 

. 12514757 
.03140785 
41  . 

270. 

.79015067 
. 79938323 

.09402992 

0.0 

.79445477 

.8 

.06276728 

.79753387 

274.5 

279. 

283.5 

288. 

292.5 

297. 

301.5 

306.0 

310.5 

315. 

319.5 

324. 

328.5 

333. 

337.5 

342. 

346.5 

351. 

355.5 

000.0 

004.5 

009.0 

013.5 

018.0 

022.5 

027.0 

031.5 

036.0 

040.5 

045.0 

049.5 

054.0 

058.5 

063.0 

067.5 

072.0 

0. 

41. 

0. 

076.5 

81  . 

85.5 

90. 

LEFT  BACK 

.8 

-.03140785.79938323 

-.06276728.79753387 

-.09402992.79445477  -.12514757.79015067  -.15607226.78462822 
-.18675629.77789594  -.21715236.76996419  -.24721360.76084521 
-.27689365.75055307  -.30614675.73910363  -.33492779.72651454 
-.36319240.71280522  -.39089699.69799681  -.41799885.68211213 
-.44445619.66517569  -.47022820.64721360  -.49527516.62825355 
-.51955844.60832477  -.54304060.58745801  -.56568543.56568543 
-.58745801.54304060  -.60832477.51955844  -.62825355.49527516 
-.64721360.47022820  -.66517569.44445619  -.68211213.41799885 
-.69799681.39089699  -.71280522.36319240  -.72651454.33492779 
-.73910363.30614675  -.75055307.27689365  -.76084521.24721360 
-.76996419.21715236  -.77789594.18675629  -.78462822.15607226 
-.79015067.12514757  -.79445477.09402992  -.79753387.06276728 
-.79938323.03140785  -.8  0.0 

41. 

270.  274.5  279.  283.5  288.  292.5 

297.  301.5  306.0  310.5  315.  319.5 
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324. 

328.5 

333. 

35 1 . 

355.5 

000.0 

018.0 

022.5 

027.0 

045.0 

049.5 

054.0 

072.0 

SITE 

076.5 

81. 

50. 

0. 

5. 

0. 

0. 

SINC 

HFLD 

pw 

0. 

1.7893 

1 . 

1  . 

8. 

1  . 

0. 

SRCS 

4. 

180. 

0. 

0. 

0. 

180. 

0. 

0. 

0. 

-90. 

-90. 

90. 

0. 

-90. 

0. 

0. 

0. 

STHE 

4. 

-90. 

90. 

1 . 

1. 

2. 

2. 

1. 

2- 

3. 

1 . 

2. 

4. 

SPRI 

1  . 

2. 

1 . 

1  . 

1  . 

1 . 

SEND 

/EOF 

1  . 

1 . 

337.5 

342. 

346.5 

004.5 

009.0 

013.5 

031.5 

036.0 

040.5 

058.5 

063.0 

067.5 

85.5 

90. 

0. 

50. 

1  . 

0. 

0. 

0 

2. 

0. 

0 

2. 

0. 

0 

2. 

0. 

0 

2. 

0. 

0 

0. 

0. 

0. 

0. 

1 . 
1 . 
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APPENDIX  E 


ORIGINAL  AND  REVISED  INITIALIZATION  FOR  VlSOL 

The  two  decks  listed  in  this  appendix  were  used  to  initialize  variables 
for  the  solver  VlSOL  source  program  AA3DS.  The  first  listing  is  from  the 
original  code  where  a  call  statement  in  AA3DS  calls  subroutine  INIT  which  in 
turn  contains  a  call  statement  to  call  up  the  BLOCK  DATA  INICOM.  The  second 
listing  is  the  result  of  solving  an  execution  error  by  imbedding  the 
initialization  in  the  source  program  AA3DS. 
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♦♦♦♦♦APPENDIX  E  -  INITIALIZATION  OF  PARAMETERS  WITH  SUBROUTINES  FROM 
THE  ORIGINAL  SOLVER  DECK  V1SOL* **♦♦♦**♦♦**»♦♦*♦**♦*♦***♦♦♦♦»♦♦♦♦*♦♦ ♦ 

♦DECK  INICOM 

BLOCK  DATA  INICOM 
C 

C  BLOCK  DATA  INITIALIZATION 
C 

♦CALL  FORMARK 
C 

♦CALL  PARAM 

C 

♦CALL  COMMARK 
C 

♦CALL  /CBPARM/ 

♦CALL  /CHRDAT/ 

♦CALL  /DIAGN/ 

♦CALL  /FILDAT/ 

♦CALL  /GRIDQ/ 

♦CALL  /INCANG/ 

♦CALL  /INCFLQ/ 

♦CALL  /ITERQ/ 

♦CALL  /MATDAT/ 

♦CALL  /NUMRCN/ 

♦CALL  /PHYSCN/ 

C 

C  INITIALIZE  FUNDAMENTAL  COMMON  BLOCK  VARIABLES 
C 

C  /CBPARM/ 

C 

C  /CHRDAT/ 

DATA  REVDAT/’ 7/31/87 7 

DATA  TITLE  / ’ ‘DEFAULT  ‘DEFAULT  ‘DEFAULT  ‘DEFAULT V 
DATA  NAME  /’ ‘DEFAULT  ‘DEFAULT  ‘DEFAULT  ‘DEFAULT 7 
DATA  CHPOL/’EE7,CHPW/’PW7 
C 

C  /DIAGN/ 

DATA  EPS/1 .E-8/, LIMIT/O/, NDOUT/6/ 

DATA  DEBUG/. FALSE 7 

C 

C  /FILDAT/ 

DATA  NIN/5/ , NOUT/6/ , NERR/6/ , NPRSL/7/ , NFCJ/9/ 

DATA  NFRCS/10/ 

DATA  NDFOUT/4/ 

DATA  NFGRN/1 1 / , IFSRCH/12/ , IFOP/13/ 

DATA  IFNJ/14/ , IFJ/ 15/ , IFGRH/ 1 6/ , IFGAM/17/ 

DATA  IFW/18/ , IFDJ/19/ 

DATA  NFMRCS/29/ 

DATA  IFQC/20/ , IFQB/21/ , IFQCI/22/ , IFQBI/23/ ,  IFQS/24/ 

DATA  NFQW/25/ ,NFQWI/26/ , IFQV/27 / 

DATA  IFQE  /30/,  IFQSC  /31/,  IFQSB  /32/ 

DATA  NFIEL/28/ 

DATA  ISEQ/1/, IRAN/2/ .IMEM/3/ 

C 

C  /GRIDQ/ 

DATA  NX  / 1 7/ ,  NY  /17/,  NZ  /17/ 
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oo  oo  no  oo 


DATA  XI  /0./,  XF  /I . / ,  YI  /0./,  YF  /I ./,  ZI  /0. / , 
C 

C  /INCANG/ 

DATA  NASW/ 1 / , NANG ( 1 ) / 1 / 

/INCFLQ/ 

DATA  AANG(l)  /0 . / ,  AANG(2)  /O . / ,  AANG(3)  /O . / 
DATA  FREQ  /3.E-8/,  EMAG  / (1 . ,0. )/ 

/ITERQ/ 

DATA  TOL  /l.E-10/ 

DATA  NITER  /30/,  ITPRT  /O/,  KITER  /20 /,  MSGLVL/3/ 

/MATDAT/ 

DATA  NMAT/2/ 

DATA  ALPHA (1) / (0 . ,0. ) / , ALPHA (2) /(O. ,0.)/ 

DATA  BETA(l) / (0. ,0. ) / , BETA (2)  / (0. ,0.)/ 

/NUMRCN/ 

DATA  CONE  / (0 . ,1.)/ 

C 

C  /PHYSCN/ 

DATA  VLIGHT  / 2.997925E8/ 

C 

END 


♦DECK  INIT 

SUBROUTINE  INIT 
C 

C  INITIALIZE  DATA 
C 

♦CALL  FORMARK 
C 

♦CALL  COMMARK 
C 

♦CALL  LOCMARK 
C 

♦CALL  CODMARK 
C 

CALL  INICOM 
C 

RETURN 

END 


ZF  /!./ 
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....♦APPENDIX  E  -  INITIALIZATION  OF  PARAMETERS  IMBEDDED  IN  SOLVER 
V1SOL  SOURCE  PROGRAM  AA3DS******»»***+*»******************»***»** 

♦DECK  AA3DS 

PROGRAM  AA3DS 
C 
C 
C 

•CALL  FORMARK 
C 

♦CALL  COMMARK 
C 

♦CALL  PARAM 
C 

♦CALL  LOCMARK 
C 

♦CALL  CODMARK 
C 

♦CALL  /CBPARM/ 

♦CALL  /CHRDAT/ 

♦CALL  /DIAGN/ 

♦CALL  /FILDAT/ 

♦CALL  /GRIDQ/ 

♦CALL  /INCANG/ 

♦CALL  /INCFLQ/ 

♦CALL  /ITERQ/ 

♦CALL  /MATDAT / 

♦CALL  /NUMRCN/ 

♦CALL  /PHYSCN/ 

♦CALL  /SCRTCH/ 

♦CALL  /STDOP/ 

♦CALL  /STENCL/ 

C 

C 

C 

CALL  CSTPRT( ’BEGIN’) 

C 

C 

C  INITIALIZE  DATA 

C 

C  BLOCK  DATA  INICOM 

C 

C  BLOCK  DATA  INITIALIZATION 

C 

C  INITIALIZE  FUNDAMENTAL  COMMON  BLOCK  VARIABLES 

C 

C  /CBPARM/ 

C 

C  /CHRDAT/ 

DATA  REVDAT/ ’7/31/87 7 

DATA  TITLE  / ’ ♦ DEFAULT  ♦ DEFAULT  * DEFAULT  ♦DEFAULT’/ 

DATA  NAME  / ’ ♦ DEFAULT  ♦DEFAULT  ♦ DEFAULT  ♦DEFAULT’/ 

DATA  CHPOL/’EE7,CHPW/’PW7 
C 

C  /DIAGN/ 

DATA  EPS/1. E-8/, LIMIT/O/, NDOUT/6/ 
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DATA  DEBUG/. FALSE./ 

C 

C  /FILDAT/ 

DATA  NIN/5/ ,  NOUT/6/ ,  NERR/6/ ,  NPRSL/7/ ,  NFCJ/9/ 

DATA  NFRCS/10/ 

DATA  NDFOUT/4/ 

DATA  NFGRN/11/ , IFSRCH/12/ , IFOP/13/ 

DATA  IFNJ/14/ , IFJ/15/ , IFGRH/16/ , IFGAM/17/ 

DATA  IFW/18/ f IFDJ/19/ 

DATA  NFMRCS/29/ 

DATA  IFQC/20/ , IFQB/21 / , IFQCI /22/ , IFQBI/23/ , IFQS/24/ 

DATA  NFQW/25/ , NFQWI/26/ , IFQV/27/ 

DATA  IFQE  /30/ ,  IFQSC  /31/,  IFQSB  /32/ 

DATA  NFIEL/28/ 

DATA  ISEQ/1/, IRAN/2/, IMEM/3/ 

C 

C  /GRIDQ/ 

DATA  NX  / 1 7/ ,  NY  /17/,  NZ  /17/ 

DATA  XI  /O. / ,  XF  /!./,  YI  /O./,  YF  /I . / ,  ZI  /O . / ,  ZF  /l./ 
C 

C  /INCANG/ 

DATA  NASW / 1 / , NANG ( 1 ) / 1 / 

C 

C  /INCFLQ/ 

DATA  AANG(l)  /O . / ,  AANG (2)  /O./,  AANG(3)  /O./ 

DATA  FREQ  /3.E-8/,  EMAG  / ( 1 . , 0 . ) / 

C 

C  /ITERQ/ 

DATA  TOL  /l.E-10/ 

DATA  NITER  /30/,  ITPRT  /Of,  KITER  /20/,  MSGLVL/3/ 

C 

C  /MATDAT/ 

DATA  NMAT/2/ 

DATA  ALPHA(l) / (0. ,0. ) / , ALPHA (2) / (0 . ,0.)/ 

DATA  BETA(l) /(O. f 0 . ) / » BETA (2) / (0 . ,0.)/ 

C 

C  /NUMRCN/ 

DATA  CONE  /(O. ,1.)/ 

C 

C  /PHYSCN/ 

DATA  VLIGHT  /2.997925E8/ 

C 

C  ♦♦♦END  OF  INITIALIZATION  DATA*+* 

C 

CALL  DATE(ADATE) 

CALL  CLOCK (ATIME) 

CALL  VHEAD(NOUT, ’V1S0L : ’ ,REVDAT) 

CALL  CSTPRT(’INIT’) 

C 

C  SET  UP  MEMORY  MANAGER 

CALL  SETRA (MAXSCR , SCRGLO) 

C 

WRITE(N0UT,2) 

2  FORMAT (5X,  ’♦♦♦JUST  BEFORE  CALL  SUBROUTINE  READING ’) 

C 

CALL  READIN(IREAD) 
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c 

WRITE(N0UT,3) IREAD 

3  FORMAT (5X,’»*»IREAD==\  15) 

C 

WRITE(N0UT,4) 

4  FORMAT (SX, ’***BACK  FROM  SUBROUTINE  READIN*** ’) 

C 

CALL  CSTPRT( ’READIN’) 

C 

IF  ( IREAD. NE.O)  THEN 
C 

WRITE(NOUT,  6) 

6  FORMAT (5X, ’*** JUST  BEFORE  CALL  SUBROUTINE  SETUP*.*’) 
C 

CALL  SETUP 
C 

WRITE(N0UT,8) 

8  FORMAT (5X , ’ *  *  *BACK  FROM  SUBROUTINE  SETUP***’) 

C 

CALL  CSTPRT(’ SETUP’) 

C 

WRITE (NOUT, 10) 

10  FORMAT (SX, ’*»♦ JUST  BEFORE  CALL  SUBROUTINE  SOLVER***’) 
C 

CALL  SOLVER 
C 

WRITE (NOUT, 12) 

12  FORMAT (5X, ’*** BACK  FROM  SUBROUTINE  SOLVER***’) 

C 

CALL  CSTPRT(’ SOLVER’) 

C 

WRITE (NOUT ,14) 

14  FORMAT (5X, ’...JUST  BEFORE  CALL  SUBROUTINE  OUTPUT***’) 
C 

CALL  OUTPUT 
C 

WRITE (NOUT, 16) 

16  FORMAT (5X, ’*** BACK  FROM  SUBROUTINE  OUTPUT***’) 

C 

CALL  CSTPRT( ’OUTPUT’) 

C 

ENDIF 

C 

STOP 

END 


57 


APPENDIX  F 


I EQUAL  DUMMY  SUBROUTINE 

This  appendix  has  the  listing  for  subroutines  CLUSS,  CLUSM  and  IEQUAL 
which  are  all  contained  in  V1S0L  source  program  AA3DS.  IEQUAL  is  called  in 
both  CLUSS  and  CLUSM  to  keep  the  memory  space  available  for  vector  array  CB 
and  KCB  intact. 
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•♦•♦♦APPENDIX  F  -  SUBROUTINES  CLUSS,  CLUSM  AND  IEQUAL  FROM  SOLVER  V1SOL 
USED  TO  PRESERVE  VECTOR  ARRAYS  IN  MEMORY* •*♦♦*♦***»♦♦♦**»*»»***♦♦*♦♦♦♦  ♦ 

♦DECK  CLUSS 

SUBROUTINE  CLUSS (N , NCB , KCB , CB , NQC , KC , NQB , KB , R , X) 

C 

*CA  FORMARK 
C 

DIMENSION  R(2.N) . X(2,N) 

DIMENSION  KCB (NCB) ,CB(NCB) ,KC(-),KB(*) 

C 

C 

*CA  COMMARK 

C 

C 

»CA  LOCMARK 
C 

COMPLEX  Y 
DIMENSION  Z (2) 

C 

C 

♦CA  CODMARK 

C 

CALL  CCOPY (N , R , 1 ,X , 1) 

KCB ( 1 ) =0 

CALL  FCHRB (NQC , KCB , NCB , KC , 1 , LC , MC , 1 ) 

DO  50  J=1,N 
NREC=KCB(1) 

IREC=J-KCB(NCB) +1 
IF (IREC.LE.NREC)  THEN 
MC=KCB(NCB-IREC) 

ELSE 

CALL  FCHRB (NQC , CB , NCB , KC , J , LC , MC , 1 ) 

END  IF 

CALL  IEQUAL (CB, KCB) 

NC=KCB(MC) 

JC=MC+NC+3 

Y=CMPLX(X(1 , J) , X(2, J) ) »CMPLX(CB(JC) ,CB(JC*1)) 

X(1 , J)=REAL(Y) 

X(2, J)=AIMAG(Y) 

IF(NC.GT.l)  THEN 
Z(1)=X(1,J) 

Z  (2)  =X  (2 ,  J) 

C  Z(l)=-Z(l) 

C  Z(2)=-Z(2) 

C  CALL  CAXPYI (NC-1 , Z,CB(JC*2) , CB(MC+4) ,X) 

C 

IK=MC*4 
M=JC+2 
CDIRS  IVDEP 

DO  25  L=1 ,NC-1 
K=KCB(IK) 

X(1,K)=X(1 ,K) -Z(l) *CB(M) +Z(2) *CB(M+1 ) 

X(2,K)=X(2,K) -Z(l) ♦CB(M+1) - Z (2) *CB(M) 

M=M+2 
IK=IK+ 1 


59 


00*00  0*0*  o  o  oooo 


25  CONTINUE 
C 

END  IF 

50  CONTINUE 
KCB(l) =0 

CALL  FCHRB (NQB , KCB , NCB , KB , N , LB , MB , 1 ) 

DO  90  1=1, N 
J=N+1 -I 

IREC= J-KCB (NCB) + 1 
IF (IREC.GE. 1)  THEN 
MB=KCB(NCB-IREC) 

ELSE 

CALL  FCHRB (NQB , CB , NCB , KB , J , LB , MB , 1 ) 

ENDIF 

CALL  IEQUAL(CB,KCB) 

NB=KCB(MB) 

JB  MB  ♦  NB  ♦  3 
IB=MB+3*NB+1 

Y=CMPLX(X(1,J) , X (2 , J) ) *CMPLX (CB(IB) ,CB(IB*1)) 
X ( 1 ,  J)=REAL(Y) 

X(2, J)=AIMAG(Y) 

IF (NB.GT. 1)  THEN 
Z(1)=X(1,J) 

Z(2)=X(2, J) 

Z(l)=-Z(l) 

Z(2)=-Z(2) 

CALL  CAXPYI(NB-1,ZICB(JB),CB(MB*3),X) 

IK=MB+3 
M=JB 
DIRS  IVDEP 

DO  75  L=1 ,NB-1 
K-KCB(IK) 

X(1 ,K) =X(1 ,K) -Z (1 ) *CB(M) +Z(2) *CB(M-1) 

X(2 ,K)=X(2 ,K) -Z (1) *CB(M+1) -Z(2) *CB(M) 

M-M+2 
IK=IK+1 
75  CONTINUE 

ENDIF 

90  CONTINUE 
RETURN 
END 

DECK  CLUSM 

SUBROUTINE  CLUSM (N , NCB , KCB , CB , NQC , KC , X , R) 

CA  FORMARK 

DIMENSION  R (2 , N) ,X(2,N) 

DIMENSION  KCB (NCB) ,CB(NCB) , KC ( * ) 


CA  COMMARK 


♦CA  LOCMARK 
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c 

DIMENSION  Z (2) 

C 

c 

*CA  CODMARK 
C 

CALL  CZERO(N,R,l) 

KCB ( 1 ) =0 

CALL  FCHRB (NQC , KCB , NCB , KC , 1 , LC , MC , 1 ) 

DO  50  J=1,N 

NREC=KCB(1) 

IREC= J-KCB (NCB) + 1 
IF(IREC.LE.NREC)  THEN 
MC=KCB(NCB-IREC) 

ELSE 

CALL  FCHRB (NQC , CB , NCB , KC , J , LC , MC , 1 ) 
ENDIF 

CALL  IEQUAL(CB,KCB) 

NC=KCB(MC) 

JC=MC+NC~3 

Z(1)=X(1,J) 

Z(2)=X(2, J) 

C  CALL  CAXPYI (NC,Z,CB(JC) ,CB(MC+3) ,R) 

C 

M-JC 
IK=MC+3 
CDIRS  IVDEP 

DO  25  L=1 ,NC 
K=KCB(IK) 

R(1,K)=R(1,K)+Z(1) *CB(M)-Z(2)+CB(M+1) 
R (2,K)=R(2 ,K) +Z(1) •CB(M+1)+Z(2)  *CB(M) 
M=M+2 
IK-IK+1 
25  CONTINUE 
C 

50  CONTINUE 
RETURN 
END 

•DECK  IEQUAL 

SUBROUTINE  IEQUAL (CB, KCB) 

RETURN 

END 
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